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ABSTRACT 

We set the ground for a theory of quantum walks on graphs- 
the generahzation of random walks on finite graphs to the 
quantum world. Such quantum walks do not converge to any 
stationary distribution, as they are unitary and reversible. 
However, by suitably relaxing the definition, we can obtain 
a measure of how fast the quantum walk spreads or how 
confined the quantum walk stays in a small neighborhood. 
We give definitions of mixing time, filling time, dispersion 
time. We show that in all these measures, the quantum walk 
on the cycle is almost quadratically faster then its classical 
correspondent. On the other hand, we give a lower bound on 
the possible speed up by quantum walks for general graphs, 
showing that quantum walks can be at most polynomially 
faster than their classical counterparts. 

1. INTRODUCTION 

Markov chains or random walks on graphs have proved to 
be a fundamental tool, with broad applications in various 
fields of mathematics, computer science and the natural sci- 
ences, such as mathematical modeling of physical systems, 
simulated annealing, and the Markov Chain Monte Carlo 
method. In the physical sciences they provide a fundamen- 
tal model for the emergence of global properties from local 
interactions. In the algorithmic context, they provide a gen- 
eral paradigm for sampling and exploring an exponentially 
large set of combinatorial structures (such as matchings in 
a graph), by using a sequence of simple, local transitions. 
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In this paper, we initiate a study of the theory of quantum 
walks on graphs — the motivation, as in the case of Markov 
chains, is to study global properties of a certain structured 
set, using repeated application of local transition rules. In 
the quantum setting, though, the local transition rule is de- 
fined to be unitary, rather than probabilistic. A classical 
Markov chain is said to be a random walk on an underlying 
graph, if the nodes of the graph are the states in S, and a 
state s has non zero probability to go to t if and only if the 
edge (s, t) exists in the graph. To define a quantum random 
walk, in addition to the Hilbert space spanned by the nodes 
of the graph, we must explicitly introduce the Hilbert space 
spanned by the outcomes of the coin that control the pro- 
cess. Thus, the quantum walk is allowed to use an auxiliary 
Hilbert space. Now, the quantum walk on a graph is nat- 
urally defined to be a unitary transformation on the tensor 
product of the Hilbert space of the graph and the auxiliary 
Hilbert space, and with the property that the probability 
amplitude (rather than the probability) is non zero only on 
edges of the graph. 

How do the basic definitions of Markov chains carry over 
to quantum walks? The most fundamental property of Markov 
chains is the fact that they converge to a stationary distribu- 
tion, independent of the initial state. However, by their very 
definition, quantum walks do not converge to any stationary 
state. This is due to the fact that unitary matrices preserve 
the norm of vectors, and hence the distance between the 
vectors describing the system at subsequent times does not 
converge to 0. One can ask whether the probability distri- 
bution induced on the nodes of the graph converges in time, 
but it turns out that it does not converge either. Yet we 
can obtain a natural notion of convergence in the quantum 
case, if we define the limiting distribution as the limit of 
the average of the probability distributions over time. This 
definition captures the amount of time the walk spends in 
each subset of the nodes, and moreover, it corresponds to 
the natural concept of sampling from the graph, since if one 
measures the state at a random time chosen from the inter- 
val {1, .., t}, the resulting distribution is exactly the average 
distribution. We show that although in general, the limiting 
distribution is a function of the initial state of the quantum 
walk, for Cayley graphs of Abelian groups it is independent 
of the initial state, and is uniform over the group elements. 

The rate of convergence, called the mixing time, is of 
crucial importance to algorithmic applications of classical 
Markov chains. Given the notion of limiting distribution 
in the quantum case, we can now talk about mixing times 
of a quantum walk. A natural definition for mixing time 
is the time it takes for the average probability distribution 
to get close to the limiting distribution. We can also talk 
about measures for how fast the quantum walk spreads or 
how long it takes the quantum walk to escape from a small 
neighborhood. We give definitions of quantum mixing time, 



sampling time, filling time, and dispersion time. How do 
the various mixing times of quantum walks compare with 
their classical counterparts? We show that the quantum 
walk on a cycle converges in time 0(n log n), giving a nearly 
quadratic speedup over the classical walk. For the cycle this 
quadratic speed up is the best possible, since the diameter of 
the graph is clearly a lower bound for the mixing time. How 
large can the quantum speed up be, for other graphs? We 
give a general lower bound on the various measures for the 
quantum mixing time, in terms of the conductance of the 
underlying graph. Our main result is that quantum random 
walks on graphs can be at most polynomially faster than 
their classical counterparts, and in fact, for bounded degree 
graphs, the gap is at most quadratic. 

It is still an open question whether quantum walks can 
be used to obtain a quadratic speed up for certain random- 
ized algorithms — such as 2-SAT. Indeed, all quantum algo- 
rithms from the last decade — including Shor's celebrated 
factorization algorithm]^ and Grover's search algorithm]^ 
— use only quantum Fourier transforms and classical com- 
putation. Is it possible to use other types of unitary transfor- 
mations to design new quantum algorithms? One constraint 
that must be met is that the unitary transformations must 
be poly-local — they must be a product of a polynomial 
number of local unitary transformations. Quantum walks 
on graphs might provide a good starting point to explore 
the effects of a sequence of local unitary transformations. 

The paper is organized as follows. We first give some back- 
ground regarding classical Markov chains and the quantum 
model. We proceed to define quantum Markov chains, and 
prove various general results about the limiting distribution. 
We then prove the speed up for the quantum walk on the 
cycle, which is followed by an upper bound on the mixing 
time for general graphs. Finally we prove the polynomial 
lower bound on the speed up for any graph, and conclude 
with a list of open questions. 

Related Work: Various researchers studied special cases 
of quantum walks on graphs. Feynman studied quantum 
walks on a line: Farhi and Gutmann[^ and Childs, Farhi 
and GutmannP studied quantum walks on various graphs 
and gave examples of graphs on which the quantum walk hits 
a particular node exponentially faster than a classical walk. 
(Note that this is a different task from the convergence to 
the stationary distribution, which we consider in this paper.) 
Ambainis, Bach, Nayak, Vishvanath and Watrous|^i studied 
various properties of the discrete-time quantum wa k on the 
line. In particular, they have calculated the the asymptotic 
behavior of the probability distributions for the walk on the 
infinite line, and shown that the probability distribution at 
time t is within a constant in total variation distance from 
the uniform distribution over an interval which is of length 
linear in t. 

2. BACKGROUND 
2.1 Random Walks 

A simple random walk on an undirected graph G{V,E), 
is described by repeated applications of a stochastic matrix 
P, where Pu,„ = J- if {u, v) is an edge in G and du the 
degree of u. If G is connected and non-bipartite, then the 
distribution of the random walk, Dt = P'-Do converges to a 
stationary distribution n which is independent of the initial 
distribution Do- For G which is d— regular, i.e. if all nodes 



have the same degree, the limiting probability distribution 
is uniform over the nodes of the graph. There are many 
definitions which capture the rate of the convergence to the 
limiting distribution. A survey can be found in [^. 

Definition 2.1. Mixing Time: 

M, = min{r| Vt > T, Do : \\Dt - tv\\ < e}, 

where here and throughout the paper, we use the total vari- 
ation distance to measure the distance between two distri- 
butions di,d2: - d2\\ = ~ d2{i)\. 

Definition 2.2. Filling Time: 

r, = min{TlVt > T,Dn,X C V : Dt{X) > (1 - e)Tv{X)}. 

Definition 2.3. Dispersion Time: 

5, = min{T|Vt >T,Do,XCV: Dt{X) < {l + e)Tv{X)}. 

The mixing time is related to the gap between the (unique) 
largest eigenvalue Ai = 1 of the stochastic matrix P, and the 
second largest eigenvalue A2. 

Theorem 2.4. Mixing time and spectral gap: / p^ 

(l-A^)log2e ^""'^ (T^K^^l^S"" + 

(1) 

The mixing time of a random walk on a graph is strongly 
related to a geometric property of the graph, the conduc- 
tance, denoted by 

Definition 2.5. Let the capacity Cx and the flow Fx 

of a subset X G G of the graph G be defined as 

Cx = ^ TTu Fx = ^ Pu,vTTu- (2) 

u£X u£X,vi^X 

where n is the stationary distribution, and pu,v is the tran- 
sition probability. Then the conductance is 



Fx 

$ — min — — 

0<|X|<|G| Cx 
Cx<l/2 



(3) 



Theorem 2.6. Conductance and spectral gap: [Jer- 
rum, Sinclairllw] 



$2 

— <(1-A2)<2$ 



(4) 



Theorems 2. A and 2.6 



together imply that the mixing time 
of a Markov chain is bounded between f2(l/$) and 0(1/$^). 

Example It is well known that for the simple random 
walk on an cycle, the mixing time is quadratic. Me = 
6i(n2-log(l/e)), and so are the filling time and the dispersion 
time. The conductance of this chain is 1/n, which gives 
a lower bound of 57(n) time steps for convergence, and an 
upper bound of 0(n^). 



2.2 Quantum Computation 

The model. Consider a finite Hilbert space Ti with an or- 
thonormal set of basis states {|s)} for s £ fi. The states 
s £ may be interpreted as the possible classical states of 
the system described by Ti. In general, the state of the sys- 
tem, |q), is a unit vector in the Hilbert space Ti, and can 
be written as \a) = Yls&n "^k), where Y.s&n l""!^ = 1- (aj 
denotes the conjugate transpose of la). {/3|a) denotes the in- 
ner product of I Of) and |/3). A quantum system can undergo 
two basic operations: unitary evolution and measurement. 

Unitary evolution : Quantum physics requires that the 
evolution of quantum states is unitary, that is the state 
\a) is mapped to U\a), where U satisfies U ■ — I, 
and [/^ denotes the transpose complex conjugate of 
U . Unitary transformations preserve norms, can be 
diagonalized with an orthonormal set of eigenvectors, 
and the corresponding eigenvalues are all of absolute 
value 1. 

Measurement : We will describe here only a measurement 
in the orthonormal basis \s). The output of the mea- 
surement of the state |q:) is an element s £ Q,, with 
probability j{s|Q)p. Moreover, the new state of the 
system after the measurement is |s). 

Combining two quantum systems : If T-La and TIb are 
the Hilbert spaces of two systems, A and B, then the 
joint system is described by the tensor product of the 
Hilbert spaces, TIa 'Si'Hb- If the basis states for Ha, 
Hb are {ja)}, {l'^)}, respectively, then the basis states 
of TLa ® Ti-B are {\a) (g) \v)}. We use the abbreviated 
notation \a,v) for the state \a) (8> This coincides 
with the interpretation by which the set of basis states 
of the combined system A,B is spanned by all possible 
classical configurations of the two classical systems A 
and B. 

Non-unitary evolution The unitary model of quantum 
computation is not the most general model possible. 
In fact, the most general quantum state is a semi def- 
inite positive trace one matrix, p, called the density 
matrix. The density matrix of |q) is |q){q:|. p evolves 
by a unitary operator U to UpU^ . In general, the evo- 
lution of the density matrix is not necessarily unitary; 
p evolves to Ep, where _B is a completely positive linear 
operator, or a super operator. Another way to think 
of non unitary evolution is by adding qubits to the 
system, applying unitary transformation on the entire 
system and then throwing the extra qubits away. For 
more details see [^, 

3. QUANTUM MARKOV CHAINS 
3.1 Definitions 

Let G{V, E) be a graph, and let Tiv be the Hilbert space 
spanned by states If) where v £ V. We denote by n, or 
\V\ the number of vertices in G. First assume that G is d- 
regular. Let TLa be an auxiliary Hilbert space of dimension 
d spanned by the states |1) through |d) (we think of this 
auxiliary Hilbert space as the "coin space"). Let C be a 
unitary transformation on TIa (which we think of as the 
"coin-tossing operator"). Label each directed edge with a 
number between 1 and d, such that for each a, the directed 



edges labeled a form a permutation. For Cayley graphs the 
labeling of a directed edge is simply the generator associated 
with the edge. Now we can define a shift operator S on TLa® 
TLv such that S|a, u) — \a,u} where u is the a-th neighbor 
of V. Note that since the edge labeling is a permutation, 
S is unitary. One step of the quantum walk is given by 
U = S • (C (8) /). We call this walk a coined quantum 
walk. 

Example: Coined Quantum Walk on the Cycle Con- 
sider the graph G which is a cycle with n nodes. This 2- 
regular graph can be viewed as the Cayley graph of the 
Abelian group Zn with the generators -1-1 (denoted by R for 
right) and —1 (denoted by L for left). The Hilbert space of 
the walk would then be (8)C". We choose the coin tossing 
operator to be the Hadamard transform, 

and the shift S is defined by 

S\R,i) = \R,i + l modn) (6) 
S|L, i) = \L,i — lmodn) 

The quantum walk is then defined to be repeated applica- 
tions of the Hadamard matrix operating on the first register, 
followed by the shift S. Note that the coin we use corre- 
sponds to a classical " unbiased" walk, in the sense that if 
measured, the walk has an equal chance of moving left or 
right. 

In our more general definition, the general quantum 
walk, we relax our restriction on the exact form of U, and 
require only that U respects the structure of the graph. 
In other words, we require that, for any v and a, the su- 
perposition U\a,v) only contains basis states \a',v') with 
v' G Q{v) U {v}, where Q{v) is the set of adjacent nodes to 
V. This means that the quantum walk only moves to neigh- 
bors of V or stays at v. More formally, let X be a set of 
vertices, and B the set of vertices which are neighbors of 
vertices in X (but that are not in X.) We denote by Px, Pb 
the probability to measure a vertex in X, B respectively. 

Claim 3.1. For any state \a) , 

Px{U\a)) < Px{\a)) + PB{\a)). 

Proof: Let |q) = |ai) -I- |q2), with |qi) being a superposi- 
tion over vertices in X U B and \a2) being a superposition 
over X — B. 

Then, PxU\a2) — because \a,v) components of |ci2) get 
mapped to components corresponding to neighbors of v and 
no vertex in X — i? is connected by an edge to a vertex 
in X. Therefore, PxU\a) = PxU\ai). Since Px is a 
projection (and can only decrease norm) and U is unitary, 
PxU\ai) < \\U\ai)\\^ = ||ai|p. Since |qi) is a superposition 
over vertices in XUB, ||ai||^ = ||Px|ai) |P -I- H-Pslai) ||^. □ 

In our most general definition, the Non Unitary Quan- 
tum Walk, we allow the quantum operation representing 
one time step of the Markov chain to be non-unitary, i.e. 
the unitary matrix U is replaced by a completely positive 
linear operator E (a super- operator) operating on the state 
of the system, represented by the density matrix p on the 
Hilbert space TL = TLc ® TLa- We say that the walk defined 
by E respects the graph G if for any density matrix p on TL 
and all subsets X of the vertices, 

Px{Ep) < Px{p)+Pb{p). 



In the rest of the paper we use the unitary definition, but 
all definitions extend in a natural way to the non-unitary 
case. 

3.2 Limiting Distribution 

We now discuss the evolution of a quantum walk as a 
function of time. Starting with an initial state |q:o), the 
state of the quantum walk at time t is \at) = l7*|ao). In 
general the limit limt^oo \at) does not exist. The reason 
being that U, as a unitary transformation, preserves the 
norm of \at) — U\at}. Consider instead the probability 
distribution on the nodes of the graph induced by \at), 

Definition 3.2. -Pt(w|ao) = I]^g^ |{a, u|at) 1^- 

We will sometimes denote this probability by P"(v). One 
might ask whether this probability distribution converges to 
a limit. However, Pt does not converge either. To see this, 
first observe that the eigenvalues of U are of the form e'^, 
and therefore after a finite number of steps, t, e*^' is arbi- 
trarily close to 1 simultaneously for all eigenvalues. Hence 
the evolution of the state is quasi periodic — the state of the 
system f/*|ao) is arbitrarily close to |ao) (and L^^^^lao) is 
arbitrarily close to |ai)) for infinitely many times t. As long 
as the probability distributions at time and 1 are difi'erent, 
Pq ^ Pi, this implies that Pt does not converge. 

Despite the fact that the actual distribution does not con- 
verge, its average over time does. We define: 

Definition 3.3. PT(t'|ao) = i ELV ■Pt("|ao) 

It turns out, as we will see soon, that for any initial state 
this quantity always has a limit as T grows to infinity, which 
we denote by 7r(w) (and sometimes write 7r(i;|ao) if we wish 
to stress its dependence on the initial vector). Intuitively, 
this quantity captures the proportion of time which the walk 
"spends" in any given node. Note that it is easy to sample 
according to this distribution Pt using the following process: 
Uniformly pick a random time t between and T —1, let the 
process evolve for t time steps and then measure to see which 
node it is at. The node will then be distributed according 
to Pt- 

We now prove a general statement about the convergence 
of Pt- The algebra used to prove this theorem will be useful 
in the rest of the paper. Let \4>i), Xj denote the eigenvectors 
and corresponding eigenvalues of U, respectively. 

Theorem 3.4. For an initial state \ao) — J2j '^jI'^j); 
lim PT(v\ao) = "S^ aiaj{a,v\(j)i){(j)j\a,v) 

Ti — >oo — ^ 

where the sum is only on pairs i,j such that Xi — Xj- 

Proof: We start by writing down the probability to measure 
the basis state \a,v) in |at), for a fixed t- 

\{a,v\at)f = I ^ a^Xl{a,v\(j>^}f 

i 

= ^ aia*{XiX*Y{a,v\(l>i){(pj\a,v) (7) 

We now take the average over time of (^), from t = to 
T—1. The only time dependent term in the above expression 



is {XiXjY . Hence, we are interested in 



T 



(8) 



We separate into two cases: One in which XiXj = 1, or 
equivalently Xi = Xj. In this case, we have that the average 
in Eq. (0) is equal to 1. In all other cases, we can write 



|iX:W*)*| = 



|1- (A,A*)' 
II- A,A* 



T\Xi Xj 



(9) 



The latter term converges to zero, therefore the contribution 
to the limiting distribution comes solely from terms with 



A, 



Xj. Thus, the limiting distribution can be derived 



from the expression in equation ^ by summing only over 
pairs which correspond to equal eigenvalues. This yields the 
desired claim, using the fact that the probability to measure 
a node w is a sum over the probabilities to measure \a,v), 
and so we can let each term converge separately. □ 

In the case in which all eigenvalues of U are distinct, the 
limiting distribution takes a very simple form. Denote by 
Pi{v) the probability to measure the node v in the eigenstate 

10,), so pi{v) = J2a \{";V\(l>i)\^. 

Corollary 3.5. If all eigenvalues ofU are distinct, then 
for an initial state |qo) ~ "^j o,j\(pj), 



lim Pt (fjoo 

— voo 



By corollary 3.5, the limiting distribution depends on the 



initial state. However, if all eigenvectors induce a uniform 
distribution over the nodes of the graph, the limiting distri- 
bution is uniform, as is easily implied by the theorem. We 
show: 

Theorem 3.6. Let U be a coined quantum walk on the 
Cayley graph of an Abelian group, such that all eigenvalues 
ofU are distinct. Then the limiting distribution n is uniform 
over the nodes of the graph, independent of the initial state 
\ao)- 

Proof: We derive an explicit expression for the eigenvec- 
tors of U, which, for a coined quantum walk, is of the form 
U — S-{C'SiI)- We note that 5* is a matrix of dimension dn, 
for n = 5 is composed of d blocks, each of dimension 
n. The a— th block corresponds to applying the a-th gen- 
erator ga on the group. We note that the characters of the 
group, Ixfe) = -^Y,„Xkiv)\v), are simultaneous eigenvec- 
tors of all the blocks. The eigenvalue associated with apply- 
ing the a-th block on \xk) is Xkida^)- Since the application 
of the coin applies an identity onTiv, a natural guess for the 
form of the eigenvectors is (J]a=i Co|a)) 'X' \Xk}- Applying 
C and then S on this vector, we find that this vector is an 
eigenvector of U if X^a'-"!'^) eigenvector of the d x d 

matrix = Ak ■ C, where Afc is a diagonal matrix, with 
Afc(a, o) = Xkida^)- Since H^, as a product of two unitary 
matrices, has d orthogonal eigenstates, the tensor products 
of these eigenstates (which depend on k) with \xk) give d 
orthonormal eigenstates for U. Running over k, this gives 
an orthonormal set of nd eigenstates for U. It is easy to see 
that the probability distribution that these eigenstates in- 
duce on the group elements is uniform, since the characters 
\Xk) are uniformly distributed over the group, and since the 



eigenstates are of the form of a tensor product, the proba- 
bihty to measure a, v in (X]a=i Ca|a)) <8) \Xk) summed over a 
is just the probability to measure v in |xfc)- This proves the 
theorem, using corollary 



i.5. □ 



We claim that for any quantum walk, if the limiting dis- 
tribution is independent of the initial node and state of the 
auxiliary space, then it must be uniform over the nodes. 

Claim 3.7. Consider a quantum walk such that for any 
initial basis state of the form \a,v), for v £ V, the limiting 
distribution over the nodes of the graph is equal to n. Then 
TT is uniform over the nodes of the graph. 

Proof: If the initial state is chosen randomly from a uni- 
form distribution over all basis states, then the limiting dis- 
tribution is equal to the average over the limiting distribu- 
tions for each initial state, but since they are all equal to 
TT, the limiting distribution for the uniform mixture is n. 
However, the density matrix which represents a complete 
mixture, i.e. a uniformly random basis state of the space 
spanned by \a, v) is preserved under unitary transformation, 
since the unitary matrix maps this space into itself. Hence 
for any time t it induces a uniform probability distribution 
over the nodes in the graph, because the initial density ma- 
trix induces this distribution. This means that the limiting 
probability distribution starting from the complete mixture 
is uniform. Combining the two facts together, we get that 
TT is uniform. □ 

3.3 Mixing Times 

We first define the analogue of the classical notion of mix- 
ing time: 

Definition 3.8. Mixing time-.The mixing time Mt, of 
a quantum Markov chain is 

= min{T| > T, \a,v) : |l7r(-|a,i;) - Pt{-\a,v)\\ < e}. 

where by the notation P{-\a, v) we mean the probability dis- 
tribution conditioned on the initial state being \a,v). This 
quantity measures the number of time steps required for the 
average distribution to be e-close to the limiting distribu- 
tion, starting from a basis state. 

We next define a closely related quantity which we call 
sampling time: 

Definition 3.9. Sampling time: T/ie SampKng time Sj, 
of a quantum Markov chain is 

S, = min{r| ^t>T,\a,v),X<ZV: 
\-n{X\a,v) ~ Pt{X\a,v)\ < en{X\a,v)}. 

This is the time it takes for the walk to approximate the 
limiting distribution point-wise. Sampling at a random time 
between and S^ — l results in a distribution which is e-close 
point- wise to the limiting distribution, justifying the term 
sampling time. In the same sense, sampling at a random 
time between and Me — 1 results in a distribution which is e- 
close to the limiting distribution in total variation distance. 

The third quantity, namely the filling time of the quantum 
Markov chain is defined as the first time at which the walk 
can claim to have visited all sets with at least (1 — e) the 
correct proportion: 



Definition 3.10. Filling time: The filling time, t^, of 
a quantum Markov chain is 

= min{r| yX CV,\a,v) 3t<T : 
Pt{X\a,v) > (l-e)7r(X|a,«)}. 

We also define the dispersion time, which is in some sense 
the opposite definition to filling time: 

Definition 3.11. Dispersion time: The dispersion time, 
^c, of a quantum Markov chain is 

Ce = min{r| VX C V, la, v) 3t < T : 
PtiX\a,v) < {l + e)7T{X\a,v)}. 

This quantity measures how fast the quantum walk escapes 
any subset of the nodes. 

Remark: We note that one could consider all the above 
definitions of mixing times with an arbitrary initial state, 
|ao), and not restrict the initial state to be a basis state of 
the form \a,v). However, the mixing time could change sig- 
nificantly. We will see in the cycle example that the mixing 
time is almost linear for initial basis states of the form \a,v), 
but it is actually quadratic for general initial states. 

The above definitions can be related one to another in 
various ways. First, it turns out that the sampling time is 
an upper bound on the mixing time, the filling time and the 
dispersion time: 

Theorem 3.12. M^,^^,n < S^. 

Proof: Fix a subset of the nodes X, and an initial state 
\a,v). Suppose at all times before Se, Pt{X\a,v) < (1 — 
e)T:{X\a,v). Then the average at time Se of the probability 
to measure X is less than (1 — e)7r(X|a, v). But by definition 
of the sampling time this is a contradiction. Hence there 
exists some time before Se at which the probability for the 
measurement outcome to be a node in X is Pt{X\a,v) > 
(1 — e)n{X\a,v), and since this is true for all X, we have 
Te < 5*5. We argue in exactly the same way to prove £,e < Se- 
The statement Me < Se follows trivially from the definition 
of total variation distance. □ 

We will later define amplified versions of these quantities, 
and find more relations between them. Let us first proceed 
to give an upper bound on the mixing time Me for the quan- 
tum walk on the cycle. 

4. QUANTUM WALK ON THE CYCLE 

In subsection |3.l| , we have defined the coined quantum 
walk on the cycle. We restrict the discussion to cycles of 
an odd number of nodes n. We first show that the limiting 
distribution for this walk is uniform. 

Theorem 4.1. The limiting distribution n for the coined 
quantum walk on the n-cycle, with n odd, and with the Ha- 
damard transform as the coin, is uniform on the nodes, in- 
dependent of the initial state |qo). 

Proof: To prov e that the limiting distribution is uniform, 
by theorem 3.6 it suffices to show that all eigenvalues of 
U are different. By the proof of theorem 3.6, the set of 



eigenvalues of U consists of all eigenvalues of the matrices: 











1 

t 

V2 











UJ 




V2 


%/2 



(10) 



where u = e~. We now show that the eigenvalues of Hfc 
are distinct. The eigenvalues of are the roots of the 
following quadratic equation: 



A"-i^sin( — )A-1 = 
n 



(11) 



The two solutions to this equation are of the form e'*'= , with 
Ok being one of the two solutions for the following equation: 



V2 



(12) 



In particular, j sin(6'fc)| < which means that the roots of 
the quadratic equation ^ are confined to two regions of the 
unit circle, Ok € [-7r/4, 7r/4] and Ok G [37r/4, -37r/4]. There 
are two solutions for equation 9k. i and 6k, 2 where 6k,i = 
TT — dk.2, so they lie in different regions, and in particular, 
they are distinct. To get equality between eigenvalues com- 
ing from different fc's, we have to have sin(2^) = sin(22^)^ 
which implies that either k — k' ot k + k' — n/2. The lat- 
ter equation has no solutions for odd n, which implies the 
theorem. □ 

Theorem 4.2. For the quantum walk on the n-cycle, with 
n odd, with the Hadamard coin, we have 



M, < 0{ 



n log n 



Proof: We prove an upper bound on the mixing time Me, 
i.e. we give an upper bound on the total variation distance 
between the average distribution Pt and the limiting dis- 
tribution TT. This is done using the following lemma which 
holds for any quantum walk. 

Lemma 4.3. Consider a general quantum walk specified 
by the unitary matrix U, and let cj)i,\i he the eigenvectors 
and corresponding eigenvalues of U , respectively. For any 
initial state \f3o) = Xli'^'ll^*)' total variation distance 
between the average probabtUty distribution and the limiting 
probability distribution satisfies 



|Pt(-|/3o)- ^(-1/30)11 



<2 E 



Proof: We recall that in the proof of lemma 3.4 we have 
already bounded the time dependent term in the average 
probability distribution. From equations (^) and we have 
that 

\Pt(v)-7v{v)\ < 



a,i,j,Xu^Xj 



T\Xi A j I 



(13) 



We now use \2ab\ < |ap + \b\'^ twice, and summing over v 
we get that \\Pt — 7r|| is at most 



E 



+ a,- 



■ \{a,v\cl>i)\^ + \{<l>,\a,v)\ 



T|A, -Aj 
(14) 



Summing first over all v and a we get the desired bound. □ 
We observe that in that lemma, the distances | Ai — Aj | are 
of crucial importance, and they need to be large for the con- 
vergence time to be small. By the proof of theorem the 



eigenvalues are distributed in two regimes (which we will call 
R and R') of the complex unit circle, 6k,i € [— 7r/4, 7r/4] = R 
and 9k.2 G [37r/4, — 37r/4] = R' . Near the boundaries of 
these regimes, i.e. for those ^'s coming from k's in the vicin- 
ity of n/4, 3n/4 modulo n, the distance between two adja- 
cent eigenvalues can be of the order of However, we 
claim that the contribution of these problematic eigenvalues 
is small, and that for the rest of the eigenvalues, the distance 
is of order 1/n. We fix < 5 < 1 (which will later be related 
to e) and define 



Rs = [0,(l-5)J]u[(l + 5)^,2^] 
R's = [(1 + 5)^(1-5)^1 



(15) 



These two regimes together cover the entire interval [0, 27r] 
except for a 25 portion of it. We refer to k such that 2it is in 
one of these regimes as "(5-good", and other fc's are "j-bad". 
We also refer to eigenvectors and eigenvalues associated with 
"(5-good" fc's as "(5-good" , and similarly for "(5-bad" . We will 
later show that if the initial state for the walk is a basis state, 
the contribution of the bad eigenvectors is small, because the 
projection of basis states on bad eigenvectors is small. But 
first, let us restrict our attention to an initial state which 
is a superposition of good eigenvectors, and consider the 
convergence to limiting distribution in this case. We first 
give a lower bound on the spacing between good eigenvalues. 

Definition 4.4. 

As — min{|Ai — Aj | s.t. i 7^ j} 

where i,j run only on S-good eigenvalues. 

Claim 4.5. For the quantum walk on the odd n cycle with 



the Hadamard coin. As > 



s/2n 



Proof: First observe that if Xi, Xj originate from the same 
k, then they lie in two different regimes R and R' , which 
means that |Ai — A^l is at least \/2- Hence, we can restrict 
our attention to eigenvalues coming from different fc's. Let 
Xi, Xj originate from fc, fc', respectively. Then using equation 
[13 we have 

lA, - A,l > 1 sm(eO - sm{ej)\ = sm( ) - sm( )| 

V 2 ri 

(16) 

We separate the proof to two cases. In the first case, fc, fc' lie 
in the same regime, Rs or R's. Recall the intermediate value 
theorem, which states that for a continuous function, for 
any x < y, there exists x < z < y such that \.f{x) — ,f{y)\ = 
\f'{z){x — y)\. Applying this theorem with f{x) — sin(a;), we 
get 

. 27rfc . ,2Tvk' 27r(fc-fc') 2tt 

sm( ) — sm( ) — cos(7)( — -)\ > cos(7) — 

n n n n 

(17) 

for some 2^ < 7 < Since fc, fc' are in the same regime, 

then -y £ Rs OT 'y £ R's, and by monotonicity of the cos 
function, we have: 

j cos(7)l > 1 cos(^ii^)l = I sin(^)l > S (18) 



where the last equahty foUows from the fact that sin(O) = 
0, sin(7r/2) = 1, and sin is convex in the regime [0, 7r/2]. If 
k, k' belong to different regimes, then we can no longer claim 
that cos(7) is large. Instead, we write 



|sin( 



2nk, 



— siri(- 



,2-Kk\ 



. . .2Tvk. . , 2-Kk\ 
— I sm( ) — sm(7r ^ 



^ I f ,27r(fc + fc') , TT 

> I cos(7 )( — ^ — TTjl > I cos(7 ) — I 



(19) 



for some 7', between — and tt— 22^1. Now, — and vr- 
lie in the same regime, and so using the same argument as 
before, the lemma follows. □ 



We can now use claim 4.5 to give a better lower bound on 



the distance between two eigenvalues. 

Claim 4.6. Let us order the eigenvalues such that < 
Arg{\i) < Arg{X2).... < Arg{X2n) < 2tt. Consider Xi and 
Xj which lie in the same regime, Rs or R'g . Then 



Xj\ > \i 

n 



Proof: Consider i,j as in the requirements of the claim. Let 
Li^j be the length of the shorter arc on the unit circle that 
connects Xj to Xi. We first claim that |Ai — A_,| > ^^Li,j 
in our regime. This is true since the ratio |Ai — Xj\/Lij 
is monotonically decreasing in 6i — 9j in the regime 6i — 
9j e [0, 7r/2], and so we can bound |Ai — Xj\/Li.j from below 
by its value on the boundary, \9i — 6j\ — n/2, which gives 



lAj — Ajl/I/ij- — 



2V2 



Hence, to bound |Ai — Aj| we give 



a lower bound on Lij. We have Lij 



-1, so it 



suffices to bound Li^+i- This is done by noticing that Li j > 
2| sinLij/2| = jAi — Aj| > A^, where the first inequality uses 
sin(a::) < x, the second equality uses simple trigonometry. 



and the last inequality uses claim 4.5 
j\Li^i+i > \i — j\As which combined with |A 
gives the claim. □ 



Hence, Li 



-A.! > 



2V2 1 



Claim 4.7. Let — '}2i'^i\4'i} such that all coefficients 
of S bad eigenvectors are zero. Then 

|,p/3_ < 2n(/n(n) + 2) 

II T II - 



Proof: Using lemma 4.3 



we write 



E I 



Tk2V2Ai 



+ 2E 



tV2 
(20) 



where in the first sum the prime indicates the fact that we 
sum over pairs i,j in the same regime Rs or R'g, such that 
Xi 7^ Xj, and in the second sum the double prime indicates 
that the sum is over pairs i,j such t hat Xi, Xj are in the 



different regimes. We have used claim 4.6 in the first sum, 
and the fact that |Ai — Aj| > \/2 for eigenvalues from different 
regimes in the second sum. 

To bound the first term, we observe that for each i there 
are at most two eigenvalues j such that |i — j| = k. We first 
sum over i,j. Then, summing over k, we use the fact that 
the sum of the first n terms of the harmonic series is less 
than ln(n) + 1. Thus, the first term is at most ^ 



For the second term we get an upper bound Using 
claim 4.5 to bound in the first term we get the desired 
claim. □ 

We now prove that the contribution of the 5-bad vectors 
is small. This follows from the following two claims. 

Claim 4.8. The projection of any basis state on the bad 
eigenvectors is of norm squared at most 25. 

Proof: The 2n dimensional Hilbert space of the quantum 
walk can be viewed as a direct sum of the two dimensional 
subspaces Lk, where Lk is the space spanned by the two 
eigenvectors originating from k. The projection of a basis 
state on Lk is of norm squared exactly The claim follows 
from the fact that there are 2Sn bad k's. □ 

Claim 4.9. Consider two initial states, jao),j/3o)- 
note by Pt ,Pt the average distributions in the quantum 
walk starting with \ao),\f3o), respectively. Then for all T, 
the total variation distance between the average distribution 
is bounded by the distance between the initial states: 



<2||lao)-lA.)|| 



Proof: Denote by |/3t) the states at time t starting with 
|Qo),|/So) as initial states. Denote by Pt' , Pf the induced 
distributions of |at), on the nodes of the graph. Clearly, 
\\Pt°'-Pt'^\\ < maxt<T WPf-PfW. Due to unitarity of the 
walk, the distance is preserved: || |at) — \Pt) \\ = \\\ao) — \Pq} II ■ 
By lemma 11 in [^, the total variation distance between the 
two probability distributions resulting from a measurement 
on two states which are e apart, is at most 2e. This proves 
the claim. □ 

Claim 4.10. Let\ao) be the initial basis state, and\Po) be 
the initial basis state projected on the 5— good eigenvectors, 
and renormalized. Then 



IP^-Tv^W <8V2S+\\P 







Proof: We write 

IIP" _ ^-11 < ||P° 



+ vr'" - n 



(21) 



The first term, by claim 4.9 is smaller than 2|||ao) — |/3o)||. 
The last term is also smaller than 2|||Qfo) — |/3o)||, since it 
is the limit of distances which are smaller than this term. 
We claim that |||ao) — I A)) II < 2^/25. This is true since by 

|ao) 



4.S we can write 



a|/3o) + l^), where \v) is a 



claim 

vector of norm at most V2S and a is larger than Vl — 26. 
Iljao) - |A))|| <\l-a\ + V25< 2V25. □ 



We can now combine claim 4.10 and claim 4.7 to finish 
the proof of the Theorem. If we now pick S = i(e/16)^, and 
T > An{ln{n) + 2)/eS, we get that 



< 8V25 + 



2n(ln(n) + 2) 



TS 



(22) 



Remark In the theory of classical Markov chains, the dis- 
tance between the first and second eigenvalues plays a cru- 
cial role in mixing time analysis. In the quantum case, we 
see that the distances between eigenvalues play a similarly 
important role; However, unlike in the classical case, since 
all eigenvalues of a unitary matrix are of absolute value 1, 
there is no special eigenvector which plays the role of the 
fixed state, and all eigenvalues play are equally important. 



5. AMPLIFICATION 

In classical Markov chains, after approaching a certain 
closeness to the limiting distribution, the distance to the 
limiting distribution starts to drop exponentially. Theorem 
4.2 gives only polynomial dependence on 1/e in the quan- 
tum case. However, one can amplify the closeness in a very 
simple way. Suppose the limiting distribution ■K{-\a,v) is in- 
dependent of the initial nod e v a nd the state a, and is equal 
to TT. (Recall that by claim 3.7 n is uniform.) In this case, 
the closeness to n can be amplified in a standard way to get 
logarithmic dependence on 1/e. This is done by running the 
walk for Me steps (i.e. for a random time between and 
Me — 1) and then measuring the node. If the measured node 
is V, we then initialize the state to be \a,v) with a random 
auxiliary state a, and start the walk again for one more stage 
of Afe steps, and so on for k times. We claim: 

Lemma 5.1. Amplification lemma Running the quan- 
tum walk for k amplification steps, each lasting A/e time 
steps, results in a distribution which is e'' close to tv. 

Proof: Define Pv.u to be the probability to measure the 
node u starting from a random initial basis state \a,v), in 
one amplification step, where a is randomly chosen from 
all basis states of the auxiliary space. The matrix P de- 
fined by these transition probabilities is a stochastic matrix. 
We claim that applying one amplification step starting from 
the uniform distribution tt one gets the distribution tt again. 
The reason is that a uniformly random state \a, v) (which in- 
duces a uniform distribution over the nodes v) is a complete 
mixture of the Hilbert space in which the walk evolves. The 
unitary transformation associated with the walk is a map 
from this space to itself, therefore, starting from a complete 
mixture of this space, the state of the system remains a 
complete mixture, i.e. vr is preserved. 

We now claim that the Li norm of any vector orthogonal 
to TT is shrunk by a factor of e by the matrix P. To prove that 
1 1 ''"-PI I ^ ^11''' II for T A- n, observe that by definition of Me, for 
any distribution a = n + t, \\a-P — it\\ < e. This means that 
for any vector r for which the sum of elements is zero, and 
each coordinate is at least — 1/n, we have \\tP\\ < e. We 
can define a basis for the subspace orthogonal to tt, which 
is composed of such vectors: These will be the vectors Vi, 
where Vi has (n — l)/n on its j"* coordinate and the rest are 
all equal to —1/n. Any r _L tt can be written as a sum of Vi 
and-Tr: r = Yliri{vi + n) = Y^^TiVi. \\tP\\ < |ril||7;iP|| < 
Ei kik = \\r\\e. 

We can now prove the claim by induction. Starting from 
a distribution a which is within S total variation distance 
from TT, we can write a = n + t where r _L tt, and ||r|| < 5. 
Then, \\aP - n\\ = ||7rP + tP - n\\ = ||rP|| < ||r||e < Se. □ 

For the cases in which the limiting distribution is inde- 
pendent of the initial state, we can now generalize our mea- 
sures of convergence to allow the possibility of amplification. 
This means that in all definitions we allow a warm start, i.e. 
we first amplify for several amplification steps, which all to- 
gether last Ta time steps, to get an initial "warm" node (the 
exact times at which one measures are chosen so that Ta is 
minimized, and the times at which the node is measured dur- 
ing Ta are referred to as the "amplification scheme"). Then 
we apply the various definitions of mixing times with the 
"warm start" node as the initial node. However, to account 
for the initial amplification stage, we add Ta to the mixing 
times. We denote those amplified versions of convergence 



with primes. Me', S!,, and so on. We have: 

Theorem 5.2. r^,^. < S', < 12E(1L^^M}1m^ 

Proof: We first prove the right inequality. Closeness to 
within e point wise will be guaranteed if the total variation 
distance is at most emmv{n{v)}. For that, by the amplifi- 
cation lemma, it suffices to apply log(l/ min7r(w))//og(l/e) 
amplification steps, each of length Me. The proof of the left 
inequality is as follows. Let Si be achieved with a certain 
amplification scheme. We then use the same amplification 
scheme for the dispersion time and the filling time, so that 
we start with the same distribution over initial nodes. Now, 
the r emaining of the proof is exactly as the proof of theorem 
|3.12| , referring only to the time interval starting at the end 
of the amplification stage. □ 

Theorem 5.3. For the quantum walk on the n-cycle, with 
n odd, with the Hadamard coin, we have 

Mi<0{nlognlog{^)) , 5^, Ce, Te' < 0(n log^ n log (i) ) 

proof: The u pper bound on S'^jT^ and ^e follows from the- 
orems and 4.2, and the fact that min{7r(?;)} = 1/n. □ 



6. GENERAL GRAPHS 

We now prove a general upper bound on Me for any quan- 
tum walk. This wi ll im plv upper bounds on the other mixing 
times by theorem |5.2| . 

Let \<j)i} be the eigenvectors of U with eigenvalues Xi. The 
upper bound will be given in terms of A, which is defined 
to be the minimal spacing between the eigenvalues. 

Theorem 6.1. Consider a general quantum walk on a 
graph G with n nodes, with an auxiliary space of dimension 
d. Then, for an initial state \a) , the total variation distance 
from its limiting distribution is 



||-PT(-|a)-7r(-|a)|| < 



7r(ln(nd/2) + 1) 
TA 



Proof: The proof follows approximately the same lines as 
the proof for the upper bound for the cycle, except that 
the complications due to throwing away part of the system 
disappear. More pre cisel y, the proof goes along the lines 
of the proof of claim 4.7, The main difference is that we 



do not have a partition of the eigenvalues i nto two regimes. 
For this reason, in the counterpart of claim [4.6| , we can have 
\Or-ej\ € [0,Tv] (instead of [0,7r/2]). ThenTlAi - Xj\/Lij 
is minimized hy \9i — 9j\ = tt (instead of 7r/2) and we get 
J A,- — \j\ > -^li — j|A. Also, the counterpart of equation 
po| has just one summation (over all eigenvalues) instead of 
two (over eigenvalues in the same regime and eigenvalues 
in different regimes). After that, we just notice that in the 
general case there are nd eigenvalues, which implies that 
k = \i — j\ runs up to nd/2. □. 

Just like we did in the cycle case, one can separate the 
eigenvectors to "good" and "bad" vectors, where the pro- 
portion of the "bad" vectors is 5, to get a better estimation 
of the mixing time. 

7. THE LOWER BOUND 

Here we are going to prove a lower bound on the vari- 
ous mixing measures of a general quantum random walk. 



In analogy to the classical case this bound will be stated in 
terms of the conductance <& of the underlying graph G (cf. 
Chapter ^). 

We will define a slightly difi'erent quantity $' first: Let 
{X,X) be a cut in the graph G (i.e., a partition of ver- 
tices into two sets). Define Bx, the boundary of X, as the 
set of vertices in X that have an edge going to X. Let 
$' = mino<lj,l<i|^l 1^. 

Theorem 7.1. The filling, dispersion, mixing and sam- 
pling times of a general quantum walk with a uniform limit- 
ing distribution are f2(l/<l>'). 

Proof: Let {X, X) be the cut that achieves the minimum 
{\Bx\/\X\). To simplify the notation, let B = Bx- 
Let Hx be the Hilbert space supported by nodes in X, i.e. 
the subspace spanned by {|a,?;)} for all v £ X and all basis 
states a of the auxiliary space. Let Tix and 7is be Hilbert 
spaces supported by nodes in X and B (defined similarly). 
We show a lower bound on filling and dispersion times by 
taking a random state |q) of form \a,v), v £ X and showing 
that the projection of (7*|q) onto Tix is small for all t < 
6(1/$')- 

For any state \(j}), let Px\4') and Pb\4'} be the projections of 

onto Hx and Hb, respectively. 
Let la) be a uniformly random basis state \a,v), v £ X. We 
bound the expected projection of \a) onto the boundary Hb 
and then use that to bound the expected projection of U*\a) 
onto TLx- (This works because the only way to go from X 
to X is through the boundary B.) 

Claim 7.2. For any t, the expected value o/ ||Pfl?7*|a) ||^ 
is at most 

Proof: \a) is a uiuformly random state in \X\d dimensions. 
Since U is unitary, (7*ja) is a uniformly random state in 
some |X|ci-dimensional subspace of H. The projection of 
this state to the |i3|d- dimensional subspace TCb is at most 
(with equality if and only ii Hb Q U^{TLx)). □ 



\B\d 
\X\d 



By applying Claim 3.1 to PxU^\oi), we get 

\\PxU''\a)\\^ < \\PxU''-'\a)f + \\PBU''-^\a)f. 



fc- 1 more times (to ||Px(7* ^1q)| 



By applying Claim 3.1 
then ||Pxi7*~^|a) II and so on), we get 

\\PxU''\a}f < ||PBt/"-»f + . . . + \\PBU\a}f + \\PB\a)\i 



By claim 7.2, the expected value of each term on the right- 
hand side (for a random |a) = \a,v}, v € X) is at most 
|-B|/|X|. Therefore, the expected value of the sum on the 
right hand side is at most fc(|_B|/|X|). For some \a) — \a,v), 
the sum J2iZo ll^s[^'|a) ||^ is at most its expectation. For 
this |q), we have 

\\PxW\a}f < J2 WPBU'Hf < J2 \\PBU'\a)f < kB. 

for all j < k. For k to be the filling time, one of ||Px[/-'|a)||'^ 



should be at least (1 — 



Then, k 



> (1-^) 



and k> (1- e)^i|j^. Since \X\ < ^\V\, \X\ > i|V| and 
k — f2(-j^) = n(l/$'). A similar argument applies to dis- 
persion time and sam pling time. The bound on mixing time 
is implied by theorem 3.12 □ 



The quantity $' (which we call boundary) is similar but 
not identical to the conductance. 

Lemma 7.3. For a graph with maximal degree d, 
<!>' < d$, 

where $ is the conductance of a simple random walk on the 
graph. 

Proof: For a simple random walk on G, the limiting dis- 
tribution is TTi, = cii,/^^d„ — dv/2\E\. Fix a cut X,X, 
and denote the conductance of this cut by "I>x. The capac- 



ity of X is Gx 



< 



d\X\ 



Let E{X : X) be the 



-2[E\ 2lEl 

set of edges going between X and X. The fiow Fx satisfies 
Fx = > 1^. Therefore, ^x = ^ > > 

IB I 

'^1^1 = This is true for any cut X, which implies the 
lemma. □ 

Therefore, r2(^) = f^(^) and theorem U implies an 
^(dk) lo'wsr bound on For constant degree d graphs, this 
lower bound is f2(^), the same as the classical lower bound 
on filling, dispersion and sampling times. Since a classical 
random walk converges in O(^) steps, this means that a 
quantum walk can be at most quadratically faster. 

Corollary 7.4. For a general quantum walk on a bounded 
degree graph, the filling, dispersion, sampling and mixing 
ttmes are at most quadratically faster than the mixing time 
of the simple classical random walk on that graph. 

For unbounded d, the factor-d gap between the two lower 
bounds (quantum and classical) is important. This gap can 
be quite large: we did not rule out the case in which the 
quantum filling time is 0(log'^ n) but d is 0(n). We suspect 
that the bound can be improved and quantum walks are at 
most quadratically faster on any graph. We can prove that 
for the special case of coined quantum walks. 

Theorem 7.5. For a coined quantum walk, the filling, 
dispersion and sampling times are S7(l/$). 

Proof: To simplify the proof, we assume that the unitary 
transformation U is of form Co S, not S oG (i.e. we first do 
the shift S and then the coin flip C). This assumption can 
be removed by replacing the starting state \a) by G~^\ot) 
and adding an extra G at the end. 

The proof is similar to Theorem We take the set X 

which achieves the conductance, $ = yu7 . Let Tlx and 
Hx be similar to the proof of Theorem |7.l| and He be the 
Hilbert space spanned by edges in the cut (i.e., the space 
spanned by \h,v), v £ X , bo v £ X). Let Px and Pc be the 
projections onto Hx and He, respectively. 

Claim 7.6. For any t, the expected value of \\PcU^\a)\\'^ 
(for a uniformly random \a) — \a,v), v £ X) is at most 

\E(X:X)\ 

\X\d 

Proof: I a) is a uniformly random state in |X|d dimensions. 
Since U is unitary, f/*|a) is a uniformly random state in 
some |X|d-dimensional subspace of H. The projection of 
this state to the \E{X : X)| -dimensional subspace He is at 



most 



\E(X:X)\ 
\X\d 



Claim 7.7. For any state \a} and any fc G N, 

\\PxU''\a)f < \\PxU''-^\a)f + \\PcU''-^\a)f. 

Proof: Define |q') = U''-'^\a). Let \a') = \a[) + lag), with 
\a'i) being a superposition over \b, v) with v £ X or bov £ X 
and IQ2) being a superposition over aU other \b,v). 
Then, 5*102) G TLx because b o v ^ X for all \b,v) that 
appear in \a2). Since the coin flip C does not change v, 
this also means that U\a'2) — C ■ S\a'2) G Ti-x- Therefore, 
PxC/|a') =^'Jft^l«'l)• 



Similarly to the proof of claim 3.1, ||PxL''|a'i)||^ < |||a'i)|| — 



\\Px\a[}f + \\Pc\a',}f = \\Px\^'}f + \\PcW)f. □ 
The rest of proof is identical to Theorem [7.l[ □ 

7.1 Lower bound for non-unitary walks 

We first deal with a special case of non-unitary walks. In 
this case, instead of U, we have a set of possible unitary 
matrices Ui, and we choose one of them randomly to apply 
at time t. The lower bounds extend trivially to this case. 

We now give a simple lower bound on the sampling time 
of a general non-unitary walk in terms of the boundary. 

Theorem 7.8. The sampling time of a general non-unitary 
quantum walk with a uniform limiting distribution are fi(^). 

Proof: Fix a cut, X,X. Suppose we start with a state 
concentrated on X. By applying the definition of a non- 
unitary quant um walk that respects the structure of the 
graph (section 3.1) several times we get 



PxE'^p < PbE*"' 



V + • • ■ + PbEp + PbP 



(23) 



For to be the sampling time, there must be some 5*6 < 
k < 25'e such that PxE''p > (1 - 3£)|§|. On the other 
hand, the sum at the right hand side of the above equation 
is exactly kPk{B\p), which, for k > Se, must satisfy 



PkiB\p) < (l + e)|B|/|G| 



This means that (1 - 3e)|^ < fc(l + e)|j5|/lG| or 
(1-3.)|X| 



2{l + e)\B\ 



(24) 



(25) 



for all X. If we pick X to be the set which achieves the 
minimum boundary |_B|/|X|, we get the desired result. □ 

We leave it as an open question to generalize the lower 
bound for non-unitary walks to other mixing measures. 

8. CONCLUDING REMARKS 

In this paper we have set up the basic definitions for quan- 
tum walks on graphs. However, the foundations of the the- 
ory of quantum walks on graphs still await discovery. We 
list here a few selected open problems. 

The first open question is for which graphs quantum speed 
up is achievable. More generally, can the 1/$ lower bound 
always be achieved quantumly? In it was shown that for 
any Markov chain, there exists a lifted version of it which 
achieves this bound, but no lifting can give better than 1/$ 
convergence. To achieve the lower bound of 1/"I> by lifting, 
one has to be able to solve the multi-commodity flow on the 
graph, a task which is in general extremely hard. There- 
fore the lifting speed-up is an existence proof, rather than 
an algorithmic one. It would therefore be very interesting 
to know whether convergence in time l/3> can be achieved 



by quantum walks for graphs other than the cycle in an 
efficiently constructible way. 

An open question is to make our two bounds tight. We 
have shown how to improve the factor of 1/d in the lower 
bound for coined quantum walks, and this needs to be gen- 
eralized to general quantum walks, or else find a counter 
example. One possible candidate is the graph consisting 
of two complete graphs connected by one edge. It is not 
clear that quantumly one cannot achieve convergence in time 
0(n) which matches the 1/d^ = l/n"I> lower bound. 

The limiting distribution for general quantum walks still 
needs to be understood. For Abelian groups, we have shown 
that coined quantum walks converge to the uniform distri- 
bution. On the other hand, we know one example in which 
a quantum walk does not converge to the same limiting dis- 
tribution as the classical simple random walk. This is a 
quantum walk on the Cayley graph of the symmetric group 
S3. Is there a simple description, perhaps via representa- 
tion theory, of the limiting distribution for quantum walks 
on Cayley graphs of non- Abelian groups? 

A very interesting question is how to use quantum walks 
in order to speed up algorithms. One way to do that is 
via speeding up the convergence time, however it is still an 
open question to give an example in which fast sampling 
cannot be done in an easy way classically. Another direc- 
tion to pursue is to find other ways of using the various curi- 
ous features of quantum walks, rather than speeding up the 
convergence time. For example, one might try to use quan- 
tum walks which converge to limiting distributions which 
are different than those of the corresponding classical walks. 
Another way might be to investigate which quantum states 
can be generated using quantum walks. Generating interest- 
ing quantum states is an important primitive for quantum 
algorithms. A well known example is the graph isomorphism 
problem which can be reduced to the problem of generating 
a certain quantum state efficiently. 
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